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Abstract 

A solution to the Einstein field equations that represents a rigidly 
rotating dust accompanied by a thin matter shell of the same type is 
found. 
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1 Introduction 



Recently, Bonnor studied axially symmetric stationary solutions of Einstein 
field equations coupled with dust and showed the reasons why a density 
gradient parallel to the axis is allowed for General Relativity but not 
for Newtonian mechanics He also considered an analytic solution 

which represents a rigidly rotating dust, and found that the solution is 
asymptotically flat in all the three spatial directions but with a total mass 
equal to zero. He attributed this to the center singularity that was believed 
to have an infinitely large negative mass, which just balances the positive 
infinite mass of the dust that fills all the spacetime. For the details, we refer 
the readers to [||]. Motivated by this particular solution, Bonnor wondered 
whether or not a non-singular rotating dust exists. In the present article, we 
obtain a solution that generates an axially symmetric rigidly rotating dust, 
which is accompanied by a rigidly rotating thin disk, namely, a singular 
hypersurface perpendicular to the axis of the rigidly rotating dust. 

It is worth mentioning that there is a paper by Georgiou 0, who 
studied rotating Einstein-Maxwell fieldsQ. Georgiou obtains exact exterior 
and matching interior stationary axially symmetric solutions for a rigidly 
rotating charged dust. His solution generates an infinitely long cylinder and 
a thin singular disk perpendicular to the axis of the cylinder. Later on 
we consider the differences and the similarities of our solution and that of 
Georgiou's. 

The present paper has also been motivated by a study by Opher, Santos 
and Wang (OSW) concerning the origin of extragalactic jets. These 
authors argued that, under certain circumstances, the spacetime given by 
what they refer to as the van Stockum metric, which is associated with a dust 
cylinder (studied by van Stockum |^ and extensively analyzed by Bonnor 

1^), can account for the coUimating effect present in extragalactic jets. 
OSW showed that this dust cylinder produces confinement for the geodesic 
motion of test particles for certain values of the radial energy and angular 
momentum. In fact, it was one of our aims to improve the OSW's model by 
looking for a spacetime that could more realistically describe a jet. 

It is worth mentioning, however, that van Stockum in fact rediscovered a 
solution which was first obtained by Lanczos 00 . Hereafter, therefore, we 
refer to the metric related to the dust cylinder as Lanczos metric. 

The extragalactic jets are ubiquitous in active galaxies, they are highly 
collimated and the matter which forms them is highly relativistic [Q. It is 

^It has been one of the referees who brought to our attention this paper by Georgiou. 
^The referees brought to our attention that Lanczos' paper preceded that by van 
Stockum. 
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worth mentioning that there is no consensus in the hterature to explain why 
they are the way they are. Many authors assume that the jets propagate 
along a direction provided by, most probably, rapidly rotating Kerr black 
holes present at the centers of active galaxies (see, e.g., Begelman, Blandford 



and Rees M and also [10, O]). This fact also suggest that putative general 



relativistic effects could be important (see, e.g., |T^, |14|, [T^, |T6|, among 
others) . 

The central engine that gives rise to the jets could be more complex than 
a simple super massive black hole, it could well occur that jets be driven by 
an axially symmetric structure present at the center of active galaxies. The 
Lanczos solution (referred to by OSW as van Stockum solution), however, 
has its weakness when applied to an actual physical situation: it represents 
an infinitely long cylinder. As a result, the spacetime is not asymptotically 
flat and has infinite mass. 

The solution here studied generates an axially symmetric rigidly rotating 
dust accompanied by a surface layer, which does not satisfy the energy 
conditions (i.e., weak, dominant and strong [|1^) in part of the hypersurface. 
We point out however, that other authors (see Refs. [jl8|, |T9|, |2^, and 
references cited therein) have investigated such structures and some of them, 
as in our case, do not satisfy any of the energy conditions. In some of 
these cases, the energy conditions may be satisfied by a suitable choice of 
parameters. 

We argue that the present study may be by itself of interest, because it 
represents a new axially symmetric dust solution, which could also motivate 
other authors to find solutions physically satisfactory. 

In section 2 a closed form of the solutions are given and the main 
properties of them are studied, while in section 3 our main conclusions and 
some discussions are presented. 

2 The Rotating Dust Metric 

Our starting point is the Lanczos metric given by 

ds^ = dt^ - 2kdtdip - Idif^ - e^{dr^ + dz^), (1) 

where 

; /• 7 2 2 2 

K = at] , T] = rt,,r , I = r — a rj , 

/^,r = ^(^,2 - ^J) , \^,z = -^V,rV,z , (2) 
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and the function ^(r, z) satisfies the Laplacian equation v^'C = 0? with 
being the Laplacian operator in Euchdean three-space. The symbol O^^; 
denotes partial derivative with respect to the argument x, and {x^} = 
{t,r,z,ip}, (/i = 0,1,2,3) are the usual axisymmetric coordinates. One 
can show that the above solutions satisfy the Einstein field equations^ 
Rfiu — gfiuRf^ = —SupUfiUu with the energy density and four- velocity of 
the dust being given, respectively, by 

P=i^2(v'. + Vl), «^ = C (3) 

The above solutions represent rigidly rotating dust. This can be seen, 
for example, by calculating the shear in this non-expanding spacetime, 
q^^ = {u^-y + Uy-^)/2, which is identically zero for the solutions given by 
Eqs.(l) and (2). However, the angular velocity of the dust, which is given by 
= {u^.y - Uy.^jl, does not vanish. 

The specific solution considered by Bonnor for example, is ^ = 
2h/ yr^^+z^ with h being a constant. As shown in [|l| this spacetime is 
free of any spacetime singularities, except for that located at the origin of 
the coordinate system, namely, r = 2; = 0. This singularity is a curvature 
singularity with an infinitely large negative mass. 

In this paper, we consider the solution with ^ = jQ{r)e^^, where JQ{r) 
denotes the zero-order Bessel function. Then, substituting it into Eqs.(2) 
and (3) we find 

k = arJi{r)e-'\ / = [l - aVi^We"'^] , 

2 

P = - — Mr)Mr)e-'^ (4) 
P = -^[Jo{r) + JKr)]e-^. (5) 

From the above equations we can see that the spacetime is singular when 
z ^ —00. To remedy this undesirable feature we can replace z by \z\ in 
Eq.(4), i.e., 

k = arJi(r)e-l^l, I = [l - a^J^{r)e-^\'\] , 

9 

P = - — Jo(r)Ji(r)e-^l^l. (6) 

Before proceeding it is worth mentioning that in a paper by Georgiou 
exact exterior and matching interior solutions are found, where the 



^In this paper we choose units such that G = 1 = c, where G is the gravitational 
constant, and c the speed of light. 
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interior solution is similar to the solution present here. As shown in 
the spacetime refers to solutions of the Einstein-Maxwell field equations for 
a rigidly rotating charged dust with vanishing Lorentz force. The solutions 
generate an infinitely long cylinder of charged dust rigidly rotating about its 
axis and a 4-current located on a singular hypersurface perpendicular to the 
axis of the cylinder at the origin of the coordinate system. 

Due to the fact that Georgiou's interior solutions present a non null 
electromagnetic 4-potencial, namely, = (0, 0, 0, A^), his equations related 
to the fi function involves the ^3 function (see Eqs.(4.2) and (4.3) in 0). 
In such a way he could set F = = 1. On the other hand, we have 
F = 1 and fi is given by Eq.(4), as a result our solutions are different. Apart 
from the constants our mass density contains in addition the term e"'^, as 
a result a stronger dependence on the z coordinate occurs as compared to 
the Georgiou's mass density. Our solution describes rigidly rotating neutral 
dust and Georgiou's solution, rigidly rotating charged dust. Consequently 
the resulting spacetimes are different. 

One can show that such resulted spacetime present here is asymptotically 
flat in the z direction, since as \z\ 00 for any particular finite value of r, 
A; — > 0, / — s> and — 1. On the other hand, the behaviour of the solution 
as r — > 00, for any particular finite value of z, shows that, for example, 
k oscillates infinitely between —00 and +00. This would indicate that the 
spacetime is not asymptotically fiat. 

On the other hand, the Kretschmann scalar is given by 

2 -2i\z\+^,) 

n ^ R^p,sR'^^' = {lQ[Jl{r)-rUr)J,{r) 

W [Jlir) + Jlir))] - 4aW2|^l [24(r) + j\{r) 
-rlJl{r)Jl{r) - ArUr)Ur) {Jl{r) + Jl{r))] 
+aVe-^'^' [Jl{r) + jf (r) + ?>Jl{r)Jl{r) {jl{r) + Jl{r))] ] 
+7eo5(^), (7) 

where TZQ{r) is a bounded function of r (see the discussions following Eq.(14) 
below), and 5{z) denotes the Dirac delta function. Using the relations 

a; ^ 0, 

(8) 

X +00, 

we can see from Eq.(7) that TZ —>■ finite, as r — > 0, and that TZ —>■ 0, as \z\ or 
r +CX3 indicating that the spacetime is asymptotically flat. Although we 
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have considered here some discussion on flatness, it is worth bearing in mind 



that some authors argue that this concept is not well defined (see,e.g., [p]] ) 
Also, from Eqs.(6) and (8) we have 



AX 



(9) 



as r — s> 0. Hence, the axis (r = 0) of the spacetime is well defined and locally 
fiat. From the above equations we can also see that, by properly choosing 
the constant a, we may have g^^p < for any r. That is, the spacetime may 
be free of any closed time-like curves. Therefore, the solution given by Eq. (6) 
represents an axially symmetric and rigidly rotating dust spacetime. 

It should be noted that the replacement of z by \z\ gives rise to a 
thin matter shell. As a matter of fact, this replacement mathematically 
is equivalent first to cut the original spacetime given by Eq.(4) into two 
parts, z > Q and z < 0, and then join the part z > Q with a copy of it 
along the hypersurface z = 0, so that the resulted spacetime has a refiection 
symmetry with respect to the surface. After this cut-paste operation, the 
spacetime is no longer analytic across the surface ^ = 0. Actually, the metric 
coefficients are continuous, but their first derivatives with respect to z are 



not. Then, according to Taub's theory [0, a thin matter shell appears 
on the hypersurface. Introducing the quantity h^^ via the relation 

where (7+^ (g^i,) are quantities defined in the region z > {z < 0), we find 
that the non-vanishing components of b^^, are given by 

bn = 622 = -2aVjo(r)Ji(r)e'^°, 

603 = feso = 2arJi(r), 633 = -4a Vj2(r), (11) 

where fiQ = fi{r,z)\z=o. Then, the surface energy-momentum tensor r^j^ is 
given by ||18[ 



(12) 

where ra^ is the normal vector to the hypersurface z = 0, given by = 6"^, 
with n = n\n^ and b = 6^. Substituting Eq.(ll) into Eq.(12), we find that 
T^i, can be written in the form 

Svrr^i, = cru^u^ + px^x^ + q{u^Xy + u^x^,), (13) 
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with 

a = —p = — a^r Jo(r) Ji(r)e~^'', 
q = aJi{r)e-^\ (14) 

where is the four-velocity of the dust restricted to the surface 2; = 0, and 

^3 



is a space-hke unit vector on the surface, given by = r5j^, and has 



the properties: X\X'^ = —1 and xxu^ = 0. The non null components of r'^^^ 
are explicit shown in the appendix, where we also obtain them using the 



alternative technique derived by Israel ^3 

Eq.(13) shows that a represents the surface energy density of the shell 
measured by observers who are comoving with the dust fluid, and p is the 
pressure in the direction perpendicular to the observers' world lines, while q 
represents the heat flow. On the other hand, Eqs.(14) show that all these 
quantities are finite for any r. One can also see that a oscillates between 
positive and negative values; when a < the energy conditions are 
violated. Note that the function 7lo{r) appearing in Eq.(7) is a combination 
of them and is finite, too. 

Note that the quantities q^u and Wfj^u defined above contain only first 
derivatives of the metric coefficients, as a result they can be written generally 
in the form 

Y,, = Y+H{z) + Y-,H{-z), 

where Y^{= {g^j,, ti'^j,}) are quantities calculated respectively in the regions 
z > and z < 0, and H{z) is the Heavside function, which is 1 for z > 0, 
1/2 for 2; = 0, and for 2; < 0. Since are all equal zero, we can see that 
g^t/ is zero even on the surface z = 0. That is, the matter shell is also rigidly 
rotating. 

On the other hand, in terms of r^i,, the Einstein field equations for the 
solution (6) takes the form 

R^lu - ^QfiuR = -S7T[pUf,Uu + T^uS{z)], (15) 

where p is given by Eq.(5) but with z being replaced by \z\. Thus, the 
solution given by Eq. (6) represents a rigidly rotating dust accompanied by a 
matter shell of the same type. 

To further study the spacetime of this solution, let us consider the total 
mass of the spacetime. Using the formula of Tolman [^], we find that the 
mass inside the cylinder r = i? is given by 

p2Tr p+00 pR 

M{R) = / {^0 + y/^[pu\o + r\5{z)]}dipdzdr 

Jo J~oo Jo 

= ^{2MR)Ji{R) + R[J^{R)-JhR)]}, (16) 
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where i? is a constant, and f^^ is the so-called gravitational energy- momentum 
pseudo-tensor [^. The combination of Eq.(8) and Eq.(16) shows that M{R) 
oscillates infinitely between — oo and +00 as i? — +00. Thus, in the present 
case the total mass of the spacetime is not well defined. Were a always 
negative it could occur that the total mass of the dust structure were zero. 
This behaviour of M{R) as i? — 00 is related to the fact that a oscillates 
between positive and negative values. To remedy this problem, one might 
cut the spacetime and smoothly match the dust structure where a > to 
vacuum spacetimes. As result one would have a hollow rigidly rotating dust 
structure with two vacuums, one inside and other outside it. It could occur 
however that this procedure generated unphysical surface layers. Due to the 
complexity of these issues, we have not yet been successful in this direction. 

3 Discussion and Conclusions 

Motivated by the article by Bonnor [|l| we have studied a particular axially 
symmetric rigidly rotating dust solution, and found that it is accompanied by 
a thin disk located on an hypersurface perpendicular to the symmetry axis. 
The undesirable feature is that the thin disk has negative energy density in 
part of the hypersurface. 

The solution we have found is in some sense similar to a solution found by 
Georgiou . This author obtains exact interior and matching exterior axially 
symmetric solutions of the Einstein-Maxwell fields equations, for rigidly 
rotating charged dust. The fact that our metric function /x be given by Eq.(4), 
instead of being /x = 1 as in [Q, turns the spacetimes different. Therefore, the 
fact that the dust be charged modify significantly the spacetime generated. 

A completely satisfactory solution for rigidly rotating dust fiuids is yet 
to be derived and deserves to be investigated. Such a solution must be 
asymptotically fiat, have finite mass, have non-singularities, and whether a 
thin shell appears it must satisfy the energy conditions. 

To see whether or not our solution allow confinement, one needs to follow a 
similar procedure as in 0. However, due to the fact that the metric functions 
now depend on both of the r and z coordinates, the study of geodesic motions 
become very complicated. We therefore leave such an issue for a future study. 
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A Alternative Calculation of r. 



Here we show that Taub's and Israel's techniques to calculate the surface 
energy- momentum tensor, t^^, are equivalent. 
We start from 



1 



Stt hm / Tf",, dz =-\im [R'i,-- K R dz 
lim / \^^^~\k^ dz (17) 



e 



[see, e.g., [J^,^). 

Due to the presence of the thin disk, i?'^ and R contain the terms 
and 91, respectively. These terms contain delta functions, and appear 
only with components of i? whose expressions contain second order partial 
derivatives with respect to z. Following 0, it is easy to show that the non 
null components of ^i^^ read 



= 2 k -I (18) 

where the functions with superscript + (— ) stand for the functions defined 
in the region z > {z < 0). 

Substituting Eq.(18) into (17) it is straightforward to show that the non 
null components of r'^^ are given by: 



87rr°o = -a^[JUr) +rJo{r)Ji{r)]e-''' 
87rr°3 = arJi{r)[l + a^J^{r)]e-^''' 

8ttt\ = Ji(r)e"^° 

r 

871t\ = aViir)-rMr)Mr)]e~f'° , (19) 

which agree with Eqs.(13), showing, therefore, that Taub's and Israel's 
techniques are equivalent. 
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Finally, it is worth mentioning that Georgiou's definition for r'^^ is a little 
bit different from Israel's. The former integrates T'^^ with respect to the 
proper distance measured perpendicularly through the thin disk. On the 
other hand, Eq.(17) is an integral of T^^^^ with respect to the z coordinate. If 
one followed Georgiou's definition, instead of having the term e~'^° in Eq.(19) 
one would have e~''°/^. 
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